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We report quadruple configuration interaction calculations within the extended 
Pariser—Parr—Pople Hamiltonian on the excited states of aggregates of polyenes, crystalline 
acenes, and covalently linked dimers of acene molecules. We determine the precise energy orderings 
and analyze the cluster wave functions in order to arrive at a comprehensive physical understanding 
of singlet fission in these diverse families of materials. Our computational approach allows us to 
retain a very large number of basis states and thereby obtain the correct relative energy orderings 
of one electron—one hole Frenkel and charge-transfer excitons versus intra- and intermolecular 
two electron—two hole triplet—triplet excited states. We show that from the energy orderings it 
is possible to understand the occurrence of singlet fission in polyene and acene crystals, as well 
as its near total absence in the covalently linked acene dimers. As in the acene crystals, singlet 
fission in the polyenes is a multichromophoric phenomenon, with the well-known 2 1 Aff playing 
no direct role. Intermolecular charge-transfer is essential for singlet fission in both acenes and 
polyenes, but because of subtle differences in the natures and orderings of the aggregate excited 
states, the mechanisms of singlet fission are slightly different in the two classes. We are thus able 
to give qualitative physical reasoning for the slower singlet fission in the polyenes, relative to that 
in crystalline pentacene. Our work also gives new insight into the complex exciton dynamics in 
tetracene crystals, which has been difficult to understand theoretically. Our large-scale many-body 
calculations provide us with the ability to understand the qualitative differences in the singlet 
fission yields and rates between different classes of 7r-conjugated materials. 


■ INTRODUCTION 

Singlet fission (SF) in organic materials, the process 
by which an optically excited spin—singlet intramolecu¬ 
lar state S* dissociates into two lower-energy spin triplet 
excitons T, has been the focus of extensive investigation 
over the past few years. m Provided specific energetic 
considerations are met in an organic donor—acceptor het¬ 
erostructure, |3j each long-lived triplet exciton gener¬ 
ated through SF in the donor molecule can in principle 
undergo charge dissociation at the donor—acceptor in¬ 
terface, thereby doubling the number of charge carriers 
compared to that obtained from the dissociation of the 
singlet exciton. SF-driven strong enhancement of per- 
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formance has been found in pentacene—C6o solar cells, 
BE] and the process can clearly have a significant im¬ 
pact on the search for organic solar cells with even higher 
quantum efficiency. BIZ] The fundamental requirement 
for efficient SF, based on energetic considerations, is 
E( S*) ^ 2 x E'(Ti), where E( S*) and E( Ti) are the en¬ 
ergies of the singlet optical exciton and the lowest triplet 
exciton, respectively. Among various classes of organic 
materials fulfilling this requirement, nua aggregates of 
tetracene and pentacene, BEIIBHIS] acene derivatives, 
[20U23] carotenoids and polyene crystals, [24H281 and 1,3- 
diphenylisobenzofuran [29] exhibit SF with T yield rang¬ 
ing from 100% to 200%. There also exist tetracene and 
pentacene derivatives which do not undergo SF. Ha eh] 
A unified theoretical understanding of SF in the various 
systems is clearly desirable for the synthesis of new ma¬ 
terials exhibiting SF and with potential use in solar cells. 
This is the goal of the current theoretical research. 

Existing theories of SF [23J [321443] have almost uni¬ 
versally focused on acenes (on pentacene in most cases), 
where it is generally agreed upon that SF is a multichro¬ 
mophoric process, and the triplets are generated from a 
triplet—triplet (hereafter 1 (TT)) state with two triplet 
excitons localized on different chromophores and with 
overall spin angular momentum of zero. The 1 (TT) state 
is also referred to as a multiexciton (hereafter ME) in 
the literature, as it is predominantly a correlated two 
electron—two hole (2e—2h) excitation from the ground 
state. The single question that has been investigated 
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most intensively in this context is whether or not charge- 
transfer (CT) excitons (hereafter CTX) mediate the S* 
to 1 (TT) conversion. [3211341 I36H391 l44l [45] Some in¬ 
vestigators have claimed that the 1 (TT) state is reached 
from S* in a single step, driven by Coulomb interactions 
consisting of two-particle intermolecular electron trans¬ 
fers. [3211341 [39] The origin of this idea lies in the es¬ 
pecially rapid SF in pentacene, where SF occurs in ^ 
80 fs. [14] According to these authors, mediation by the 
CTX, being necessarily a two-step process (S* CTX 
—>> 1 (TT)), would make 1 (TT) generation and the SF fol¬ 
lowing it much slower. This viewpoint has been contra¬ 
dicted by the proponents of the CTX-mediated mecha¬ 
nism, [lk [36h38l l44l 145] who believe that the only require¬ 
ment for the two-step process to be ultrafast is that S* 
and 1 (TT) are of nearly the same energy. Additional jus¬ 
tification of the CTX-mediated process comes from the 
matrix element calculations of the Coulomb interaction 
involving two-particle hops, which have found this matrix 
element too small to explain femtosecond ME generation. 

Theories emphasizing CT also differ in details them¬ 
selves. A number of them find the CTX significantly 
higher in energy than S*, [Ml EH 3S1 39] in which case a 
superexchange like mechanism of 1 (TT) generation from 
S* has been suggested. Yet other mechanisms have also 
been proposed. Time-resolved experiments have been in¬ 
terpreted as showing direct photoexcitation of a coherent 
superposition of S* and X (TT). [12] Yamagata [46] and 
Beljonne et ah, [38] based on their fitting of the Davy¬ 
dov splittings in the absorption spectra of the acenes, 
claim that optical excitation here is to a state that is a 
quantum mechanical superposition of S* and CTX and 
that transition from this optically accessible state to the 
1 (TT) occurs in a virtual one-step process. These authors 
have however also found that for intermolecular sepa¬ 
ration smaller than the equilibrium distance (i.e., very 
large intermolecular hopping) direct optical excitation to 
a state that is a superposition of S*, CTX, and X (TT) is 
possible. Zeng et al. [42] arrive at the same conclusion, 
but once again, only for artificially large intermolecular 
couplings. Finally, it has been suggested that SF occurs 
directly from S*, without going through any real or vir¬ 
tual intermediate state. m 

The above disagreements are largely due to the absence 
of understanding of the precise natures and relative en¬ 
ergies of the electronic states relevant to SF. Equally im¬ 
portantly, there exist several additional questions beyond 
those that have been investigated so far. First, how do 
morphology and environment influence SF? It has been 
widely observed experimentally that structural arrange¬ 
ment of the chromophores often determines whether or 
not SF efficiency is high, [l] [2j [28l [31] Second, what 
is the mechanism of SF in carotenoid aggregates, m v 
[28] in which the well-known intramolecular 2 1 A“ state, 
which is a quantum-entangled single-molecule state of 


two triplets, [471149] occurs below the optical 1 X B+ exci- 
ton? [50] [51] (Here the superscript 1 indicates spin sin¬ 
glet nature of the eigenstate while the “plus” and “mi¬ 
nus” superscripts refer to the charge-conjugation sym¬ 
metry). Is SF here monomolecular, with the triplet ex¬ 
citons originating from the dissociation of the 2 X A“, or 
is the mechanism of triplet generation the same as in 
acenes? In either case, why is SF in carotenoids signif¬ 
icantly slower [24H271 [52] than in crystalline pentacene? 
Third, as has been noted by several research groups, al¬ 
though E( S*) < 2 x E(Ti) in tetracene, [53] a simple ac¬ 
tivated mechanism for SF appears to fail. 16 19 What 
precisely is behind this complex exciton dynamics? Fi¬ 
nally, beyond the simple acenes and their derivatives, SF 
has also been experimentally studied in the so-called co¬ 
valently linked acene dimers. [SOI IS] However, unlike 
acene crystals or even disordered films, m little to no 
SF occurs in the dimer molecules. mm There have 
been few, if any, theoretical discussions of SF in cova¬ 
lently linked acene dimers; ia m whether or not the 
low yield of triplets in them can be understood within a 
comprehensive theory of SF remains an open question. 

Our goal in the present work is to arrive at a broad 
theoretical framework for SF that can give qualitative 
and, wherever possible, semiquant it at ive answers to the 
above questions. Unlike existing quantum chemical ap¬ 
proaches to SF, we make no attempt to calculate tran¬ 
sition rates between excited states. This is because esti¬ 
mates of transition rates depend heavily on the calculated 
relative energies of initial and final states, which include 
the one electron-one hole (le—lh) optical exciton and the 
CTX, as well as the 2e—2h 1 (TT) state and the 2 X A“. 
Treating le—lh and 2e—2h excitations on equal footing 
remains difficult for molecules with more than 8—10 elec¬ 
trons within first principles approaches. Getting the cor¬ 
rect energy ordering of 2e—2h excitations in general re¬ 
quires configuration interaction (Cl) with up to quadru¬ 
ple excitations (quadruple Cl, hereafter QCI) from the 
ground state. [54H56] Over and above in the present 
case such calculations have to be performed for aggre¬ 
gates (at least dimers) of large molecules and lie out¬ 
side the scope of standard quantum chemical techniques. 
Theoretical approaches to SF therefore mostly involve 
calculations based on a few excitations across a limited 
number of frontier molecular orbitals (MOs) of pairs of 
molecules. We will show that such approximations can 
lead to errors in the quantum mechanical descriptions 
of the multichromophore eigenstates relevant in SF. For 
our modeling of chromophore clusters we therefore take 
the opposite approach here and employ the semiempiri- 
cal Pariser—Parr—Pople (PPP) 7r-electron Hamiltonian. 
[571 EH] Our approach allows us to perform Cl calcu¬ 
lations retaining up to quadruple excitations across a 
very large number of MOs (including the complete set 
in some cases), thereby ensuring accurate cluster excited 
state energy orderings. If one now makes the reasonable 
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assumptions that (i) efficient transitions can occur be¬ 
tween energetically proximate states and (ii) for SF to 
occur the 1 (TT) should be either the lowest excited state 
of the cluster or must be close in energy to an optical 
state, then a physical understanding of SF can be ob¬ 
tained from the excited state energy spectrum alone. 

The PPP Hamiltonian has successfully addressed cor¬ 
relation effects in a variety of single chromophore systems 
[59l(6T] and in recent years has also reproduced experi¬ 
mental features in multichromophore systems. [3} (62l 63 
In the present case, we will for the first time give a precise 
explanation of why SF in polyenes is multichromophoric, 
and the 2 X A“ plays no role. A natural explanation of 
the much slower SF rate in polyenes compared to pen- 
tacene also emerges from our work. New insights on SF 
and exciton dynamics in tetracene, which have remained 
considerably less understood than in pentacene, are ob¬ 
tained. Finally, we give a clear explanation of the very 
low SF yield in covalently linked acene dimers. It is only 
through high-order Cl studies of the PPP model that 
such a unified picture is obtained. 


■ THEORETICAL MODEL, COMPUTATIONAL 
METHODS AND PARAMETERIZATION 

Model Hamiltonian. We report QCI calculations for 
clusters of chromophores based on the modified PPP 
Hamiltonian consisting of intra and interchiomophoie 
terms 
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where c^ ia creates a 7r-electron of spin a on carbon (C) 

atom i located in chromophore /q h^ icr is the 

number of electrons of spin a on atom i within chro¬ 
mophore /i, and n^i = The intramolecular 

one-electron hopping integrals are between nearest- 
neighbor C atoms i and j. We have, however, chosen 
long-range intermolecular hopping integrals tjj between 
C atoms i and j (see below). U is the Hubbard repul¬ 
sion between two electrons with opposite spins occupying 


the same atomic p z orbital; Vij is the long-range intersite 
Coulomb interaction between two electrons on a single 
chromophore; and is the corresponding interaction 
for different chromophores. 

Computational Approach. We use a basis of local¬ 
ized Hartree—Fock (HF) MOs of the individual molec¬ 
ular units in our calculations. [3] 1621(64] This allows a 
clear distinction between excitations localized on a sin¬ 
gle molecule versus those that are delocalized over mul¬ 
tiple molecules, as well as quantum mechanical super¬ 
position of these extreme configurations. For the very 
largest systems we are forced to ignore some of the out¬ 
ermost HF MOs in our QCI calculations. The number of 
MOs we have retained in all cases (cf. Supporting Infor¬ 
mation) is however significantly larger than what are re¬ 
tained within existing calculations. [33H361 l39l (42l(45l [65] 
For precise classification of all excited states proximate 
to the optical exciton in a multichromophore system, the 
QCI wave function of the nth excited state is expanded 
in terms of single, double, triple, and quadruple e—h ex¬ 
citation configurations 

l*S> = - lh >, + - 2 k>, 

i i 

+ X 4 n,3) l 3e - 3h > fe + E C ,4) l 4e - 4 A ( 4 ) 

k l 

where within each class the |^|’s are sorted in a de¬ 
scending order. For each eigenstate we now calculate 
the wave function character density p = X^|>o l 161 2 ? 
corresponding to each class of excitations i (thus the cal¬ 
culated p for a given eigenstate allows us to classify it 
as predominantly le—lh, predominantly 2e—2h, ... etc.). 
In what follows we ignore configurations with |^| < 0.1 
and also 3e—3h and 4e—4h excitations (3e—3h and 4e—4h 
excitations in any event have |^| < 0.1). 

On the basis of such detailed analyses of the wave func¬ 
tions we are able to further classify f62h64] eigenstates 
dominated by le—lh excitations as (i) Frenkel exciton 
(hereafter EXC), if the excited electron and hole are lo¬ 
calized on a monomer, (ii) polaron pair (hereafter P + P - ) 
if the electron and hole are completely charge separated 
and lie on different molecules, and (iii) CTX, if the wave 
function is a superposition of EXC and P + P - . Sim¬ 
ilarly, eigenstates dominated by 2e—2h excitations can 
have both excitations on the single monomer or on dif¬ 
ferent monomers. In the energy region of interest (near 
and below the optical exciton), 2e—2h excitations local¬ 
ized on a monomer are the predominant constituents of 
the 2 1 A~ as well as of optically dark 1 B“ states (with 
additional contributions coming from configurations that 
are EXC). [471(49] These excited states are quantum- 
entangled triplet—triplet states, where both triplet ex¬ 
citations occupy the same molecule. In the language of 
valence bond theory, these spin singlet excited states are 
covalent , ga m] as opposed to the ionic EXC states 


4 


Scheme 1: Nomenclature for the le—lh and 2e—2h Excitations that Dominate Eigenstates Energetically Proximate to 
the Optical Exciton in a Two-Chromophore (Cl and C2) System: (a) EXC, (b) P+P - (c) COV and (d) (TT) n a 
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“Excitations across the HOMO and LUMO only are shown, while actual excited states may have excitations involving much lower 
bonding and higher antibonding MOs. Covalent excited states such as the 2 1 AL^ have contributions from configurations of the type (c) 
as well as EXC configurations. 


(2e—2h states can also be singlet—singlet; in this case 
they are doubly ionic, with two pairs of C + and C _ ions; 
these have energies nearly twice that of the optical exci¬ 
ton m and are not of interest here). We will label such 
monomer 2e—2h l A.~ and 1 B~ excitations collectively as 
COV. Finally, there exist also eigenstates dominated by 
2e—2h excitations in which two le—lh excitations occur 
on different monomers. These are overall spin singlet 
excitations consisting of two triplets on two monomers, 
which we will label as (TT) n , where n is a measure of 
the separation between the triplet excitations (n = 1 
corresponds to nearest-neighbor triplets, n = 2 as next 
nearest neighbor, ..., etc.). In Scheme Ji— d we have given 
the schematics of EXC, P + P - , COV, and (TT) n excita¬ 
tions, respectively, for a dimer of two chromophores Cl 
and C2, where, however, we have included only the high¬ 
est occupied and the lowest unoccupied MOs (HOMO 
and LUMO, respectively) in describing them. The ac¬ 
tual wave functions can have contributions from excita¬ 
tions across much larger one-electron gaps. Note that 
the (TT) n configuration is a normalized superposition of 
three fundamental configurations, and the formal expres¬ 
sion in the limited HOMO—LUMO basis of Scheme [] is 




l lLa a lH-cr a 2L-cr a 2H(j 


o( a lZ4 a l Hi ~ a lLt a lHt)( a 2Li a 2Hi ~ a 2L^ a 2H^) \&) (5) 


where a\ Hcr (at Lcr ) creates an electron in the HOMO 
(LUMO) of molecule i = 1,2 and | G) is the HF ground 
state. There are three distinct terms (since a =t, i) 
which correspond to ( SS z ) = (+1,-1); (—1,+1); and 
(0, 0) excitations on the individual units, where S l z is the 
z component of the spin on molecule i. Using the same 


limited HOMO—LUMO basis, the other three classes of 
excitations in Scheme D are written as 

EXC = - ^ ( a iL(j a lHcT =*= a 2Lcj a 2Hcr) \G) 

(6a) 

P + P~ = \Y, ± a\ La a 2Ha )\G) 

(6b) 

COV = (aIz,t a lZ4 a 1 *U a 1 fft ^ a 2Lt a 2Ll a 2Hl a 2Ht)\C) 

(6c) 


Only the “plus” superposition in eq[6f is dipole allowed; 
the “minus” superposition is an optically “dark” exciton. 
CTX states are superpositions of EXC and P + P - with 
varying relative weights. 

For further clarifications of our labeling scheme we 
note the following. First, multiple-chromophore systems 
are expected to have multiple excitonic states that are 
dip ole-allowed, which is the origin of the experimentally 
observed Davydov splitting in these systems. Second, 
however, not all EXC states are optically allowed. Par¬ 
ticularly when the two molecules in a dimer cluster are 
related by any symmetry operation, the odd superposi¬ 
tion of the excitations on two molecules that are other¬ 
wise equivalent will be optically dark. 

Parameters. Our parametrizations of Hi ntra are as 
follows. The intramolecular hopping integrals are —2.4 
eV for the phenyl C—C bonds unless otherwise stated and 
—2.2 (—2.6) eV for the single (double) bonds. These val¬ 
ues are accepted as standard within correlated-electron 
models l47Tf49] (noninteracting Hiickel models often use 
larger hopping integrals —2.7 to —3.0 eV in order to 
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fit experimentally observed optical absorptions; such 
large values are neither appropriate nor necessary when 
Coulomb interactions are nonzero). The intramolecu¬ 
lar intersite interactions are obtained from parametriza- 
tion of the PPP model obtained previously, [60, Vij = 

U/k^Ji + 0.6117where Rij is the distance in A be¬ 
tween C atoms i and j and k is an effective dielectric 
constant. For the linear polyene calculations reported 
below, we use U = 8.0 eV and k = 2. These parame¬ 
ters have been used extensively in the past, mm in¬ 
cluding for linear polyenes, [56 in all cases with excel¬ 
lent fits to experimentally determined singlet and triplet 
excitation energies. The Coulomb interaction parame¬ 
ters that we chose for our calculations on the acenes are 
slightly smaller, in order to reproduce the experimentally 
obtained energies of the spin triplet excitons (see below). 

For Hi n ter, we have chosen the same functional form 
for the intermolecular Coulomb interactions as in 
Hintra • The value of the intermolecular screening param¬ 
eter k, 1 - is uncertain; we have taken this to be the same as 
the intramolecular screening parameter k for simplicity. 
The intermolecular hopping integrals between C atoms 
i and j are obtained by adjusting /3 in the expression 

[62j[66] 

tij = P e XP [(hmin ~ fi»j)/^]( ni D^ )( n2 n^ J ) (7) 

-ttij -H'ij 

with h m in being the distance between the closest two 
C atoms located on two separate chromophores, (3 the 
hopping integral between those two C atoms, and S = 
0.045 nm. Unit vectors hi and 112 are normal to the 
surfaces of the chromophores 1 and 2 on atoms i and 
j, respectively. For our calculations on polyenes, we set 
/3 — —0.2 eV, which has been previously used successfully 
for multiple-walled carbon nanotubes. |66j For acenes, we 
have performed calculations for the range of 0 ^ \/3\ ^ 0.2 
eV (cf. Supporting Information for details) to reproduce 
known experimental results. As mentioned above, and as 
we shall see from our calculations below, the parameters 
we have chosen can be further justified a posteriori , based 
on our ability to explain SF and associated trends in all 
the compounds we consider. 

■ RESULTS AND ANALYSIS 

Linear Polyenes. In carotenoids, the linear polyene 
backbone of alternant single and double bonds accounts 
for many of the physical properties related to their ex¬ 
cited states. m In what follows, the benzene rings at 
both ends of diphenylpolyene molecules have been omit¬ 
ted to reduce the number of C atoms in our molecular 
clusters to what is numerically tractable. QCI calcula¬ 
tions were performed for clusters of up to four octate- 
traene and three decapentaene molecules. We show the 


results only for octatetraene, as the results for decapen¬ 
taene were similar in all cases. We investigated three 
different cluster geometries: (a) eclipsed, in which the 
polyene molecules are symmetrically stacked on top of 
each other (Figure [lji) , (b) slip-stacked, with molecules 
from the second layer onward shifted rigidly along the 
molecular axis (Figure EF and (c) monoclinic (Fig¬ 
ure Eh The eclipsed and slip-stacked geometries are 
unrealistic; our reason for adopting them is because they 
allow us to visualize the effects of aggregation in a sys¬ 
tematic manner, by adding monomers one at a time. 
As seen below, calculations with these structures also 
help in explaining why efficient SF may not occur even 
when energetic considerations are satisfied. Experimen¬ 
tally, carotenoid crystals have either the monoclinic or 
the orthorhombic structure, (68j [69] and the former is 
associated with efficient SF. [28] 

Eclipsed Geometry. Our computational results for this 
case are summarized in Tableland Figure [2] Table [] lists 
the energies of different kinds of excited states, excluding 
the P + P - and CTX excitations. The excited states are 
identified from comparisons of the energies of the single¬ 
molecule and multiple-molecule eigenstates. These iden¬ 
tifications are further confirmed from the analyses of the 
multiple-molecule wave functions, as done in Figure [2j In 
every case in Figure [2] p gives the normalized weight of 
particular excitation components within the wave func¬ 
tion of excited states indicated in the figure (thus, for 
example, the contribution by the (TT) n excitation com¬ 
ponents to the lowest excited state of two molecules of 
octatetraene in Figure [ 2^1 is 0.72. In reality the (TT) n 
contribution would have been even larger had we not 
ignored 3e—3h and 4e—4h excitations in our computa¬ 
tion of p. Similarly the next two eigenstates in the two- 
chromophore cluster in Figure are derived from the 
single molecule 2 X A“ state and so on). States 4, 5, and 
6 in Figure[2^ are seen to be of mixed character. Thus, 
state 4 is CTX, with P + P - contribution of ^ 0.5 and 
EXC contribution of ~ 0.25. Similarly state 6 is ^ 0.55 
(TT) n and ~ 0.05 P + P - and so on. 

A number of features and trends are noteworthy re¬ 
garding the results presented in Table [] and Figure [2j 
First, there is little to no mixing between eigenstates that 
are predominantly le—lh and those that are predomi¬ 
nantly 2e—2h when the number of molecules is even. The 
absence of coherent superposition of EXC and (TT) n , 
seen with even numbers of chromophores, persists for 
other geometries and other chromophores (see below). 
Of course such a superposition can be obtained for tjj 
much larger than assumed in our work, but as we will 
see our parametrized tjj with /3 = —0.2 eV already gives 
Davydov splittings larger than that observed experimen¬ 
tally and hence gives the upper limit to realistic tjj . The 
optical excitons are easily identified in all cases from their 
nonzero transition dipole couplings with the multichro- 
mophoric ground state, as is indicated in Table [] As 
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(c) 3D monoclinic lattice-structure 
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p = 96°, 57’±10’ 
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Figure 1: Polyene clusters considered in this work, (a) Eclipsed, in which the molecules lie in parallel planes, separated by 3.5 A; 
each atom lies directly on top of another in the next plane, when viewed vertically. The nomenclature for the intermolecular 
triplet—triplet states (TT) n is defined in the figure, (b) Slip-stacked, generated by sliding each molecule in (a) from the second 
molecule onward by half the molecular length along the direction of the molecular axis. The (TT) n states are defined as in (a). 
An alternate zigzag slip-stacked geometry, in which each atom of the second neighbor molecule, but not the nearest neighbor, 
lies directly on top of the corresponding atom of the original molecule, gives very similar results as for this case, (c) The four 
molecules that constitute the nearest and next nearest neighbors of one another within the monoclinic crystal structure of 
carotenoid crystals. The molecules occupy the vertices of the gray parallelogram. The separations between the molecules, as 
well as the different (TT) n states, are indicated in the figure below. Lattice constants and angle /3 have been taken from refsl68l 

and EH 


Table 1: Calculated Energies in eV of the Optical Exciton, Ti, Lowest (TT) n and COV States of Clusters of Octatet- 
raene Molecules® 


chromophore number 

optical exciton ( p ) 

Ti 

(TT)n 

2 J A- 

I'B- 

2 (eclipsed) 

4.94 (2.35) 

1.79 

3.65 

3.77 

4.67 



1.83 


3.81 

4.71 

3 (eclipsed) 

4.72 (1.88) 

1.69 

3.58 

4.00 

5.18 


4.74 (2.46) 

1.75 


4.07 

5.26 



1.78 


4.15 

5.28 

4 (eclipsed) 

4.41 (2.91) 

1.62 

3.42 

3.98 

5.84 


4.47 (2.39) 

1.67 


4.04 

5.87 



1.72 


4.11 

6.00 



1.74 


4.20 

6.02 

4 (slip-stacked) 

4.19 (3.65) 

1.70 

3.57 

4.16 

5.94 


4.20 (1.60) 

1.71 


4.19 

5.97 



1.73 


4.22 

6.01 



1.73 



6.07 

4 (monoclinic) 

4.32 (3.88) 

1.73 

3.63 

4.18 

5.92 



1.73 


4.18 

5.97 



1.73 


4.18 

6.01 



1.73 


4.26 

6.04 


®The single-molecule energies are E( Ti) = 1.65 eV, E(2 1 A^) = 3.42 eV, ^(l 1 BL) = 4.28 eV and E( 1 1 B^) = 4.50 eV. The numbers in 
parentheses in the first column are the transition dipole couplings between the optically allowed states and the ground state, in eA. In 
the slip-stacked geometry there are no pure g or u states and the 2 1 A^‘-derived states are weakly dipole coupled to the ground state. 
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seen in Figure [2] among eigenstates with significant EXC 
character, the highest energy states are allowed optically, 
indicating the H-aggregate character of the eclipsed ge¬ 
ometry. The energy of the optical exciton decreases with 
increasing number of molecules, a well-known effect that 
is in agreement with the experimentally observed red- 
shift of the optical state in the solid. Our calculations 
also reproduce qualitatively the Davydov splitting be¬ 
tween optical excitons expected in clusters. This is the 
energy difference between S* and S 2 in the figures. 

Our criterion here for determining that a given mor¬ 
phology is conducive to SF is that the lowest excited state 
necessarily be (TT) n (see however below for discussion of 
the tetracene excitation spectrum). Of particular interest 
in the present context therefore are the relative energies 
of the (TT) n , COV, and the CTX eigenstates. We have 
found that the multichromophoric (TT) n states not only 
include the lowest triplets Ti but also can be composed 
of the next higher triplet T 2 (this is why the number of 
(TT)i states in Figure]^ exceeds 2). The COV states are 
in all cases higher in energy than the (TT) n . These are 
derived from both the 2 1 A~ and the 1 1 B“ in Figure 2^i 
and from the 2 1 A~ alone in Figure |2 Jd, c. As seen in Ta¬ 
ble [J interchain interactions increase the energies of the 
COV states 2 X A“ and 1 X B“, while the shift in the ener¬ 
gies of Ti as well as (TT) n is very weak. The increase in 
energies of the COV states, relative to the lowest (TT) n , 
has obvious significance for SF. We have therefore exam¬ 
ined this very carefully, by performing QCI calculations 
for clusters of hexatrienes, for which for up to three chro- 
mophores all MOs could be retained in the calculations. 
Results identical to those shown in Table]] were found (cf. 
Supporting Information). Furthermore, with increasing 
number of chromophores, we have ascertained that the 
(TT) n states are superpositions of n = 1, 2 and n = 3, in¬ 
dicating the tendency to give a (TT) n band in which the 
two triplets are not bound. While all of these suggest 
efficient SF in this geometry, as indicated in Figure [2j 
the optically dark CTX states progressively shift down¬ 
ward in energy with an increase in the number of chro¬ 
mophores and constitute the lowest-energy state for the 
four monomer case. The hypothetical eclipsed geometry 
is thus not conducive for SF, as the lowest excited state 
here is an excimer (the CTX is the lowest excited state of 
three eclipsed decapentaene chromophores in our calcula¬ 
tions, indicating the role of chromophore size in addition 
to the geometry). 

Slip-Stacked Geometry. We show the results for the 
slip-stacked geometry for four octatetraene molecules 
in Figure If Even for small inter molecular hopping, 
the intramolecular center of inversion is lost, and eigen¬ 
states are no longer pure g or u. A pure one-photon 
forbidden 2 X A~ state is thus replaced with a weakly 
dipole-allowed 2 1 A“-derived state. All optically accessi¬ 
ble states are strong superpositions of le—lh and 2e—2h 
configurations. More importantly, the lowest eigenstates 


are now pure 1 (TT) (the centers of inversion between 
pairs of molecules continue to persist for rigid transla¬ 
tional shifts), and the excimer states occur above these. 
The slip-stacked geometry is thus conducive to SF. While 
our conclusion here is the same as in ref [2] the physical 
reason behind our result is very different. The previ¬ 
ous results [2] were based on symmetry characteristics 
of the HOMO and LUMO of pentacene. Our result is 
a bandwidth-induced effect that is unrelated to sym¬ 
metries of the molecules or the frontier MOs; specifi¬ 
cally, large intermolecular hopping, as would occur for 
the eclipsed geometry, favors the excimer over 1 (TT). 
Note that this conclusion emerged here only after calcu¬ 
lations for four monomers were done and calculations 
based on two monomers would have predicted the incor¬ 
rect result. The smaller overall intermolecular hopping 
between MOs, as would occur for the slip-stacked geom¬ 
etry, reverses the energy ordering. 

Monoclinic Geometry. Our calculations for the mon¬ 
oclinic geometry were done for four octatetraene chro¬ 
mophores. The wave function analysis corresponding to 
the states below the optical gap is presented in Figure]^. 
Although we used the same (3 for the interchromophore 
hopping integral (cf. eq [7| as for the eclipsed geome¬ 
try, the consequences of CT here are weaker, because of 
the larger intermolecular separation as well as the de¬ 
pendence of tjj in eq[7]on the relative orientations of the 
chromophores (cf. Supporting Information for details). 
The optical state, as well as several dark states below it, 
is now almost purely EXC in character. H-aggregate be¬ 
havior is seen again, with the highest among these being 
the optically allowed excitation from the ground state. 
The results in Figure are different from those for the 
eclipsed geometry in two major ways. First, the (TT) n , 
CTX, EXC, and COV states now occur in bunches, un¬ 
like in the eclipsed case, where the energetic locations of 
states with the same character were more dispersed. Sec¬ 
ond, the lowest six states are now all (TT) n in character 
and are separated from the COV states (all derived from 
the 2 1 A~ in this case) by P + P - excitations. The (TT) n 
states are again at energies slightly higher than 2 x E{ Ti) 
(which implies exoergic separation into triplets) and also 
show signatures of band formation (indicating that the 
triplets are not bound). Efficient SF, in agreement with 
experimental observation [28], is expected here. 

The results of Figure If are in strong support of (i) 
the multichromophoric nature of SF in the polyene and 
carotenoid molecular crystals with delocalization of the 
triplets to distant chromophores due to aggregation and 
(ii) mediation of SF by the intermediate P + P - or CTX 
states. The mechanism of SF that one envisages from 
Figure [3 is rapid nonradiative relaxation from the op¬ 
tical exciton to the lowest (TT) n state, via the interme¬ 
diate states. Whether or not the covalent 2 1 A~ states 
take part in the SF is difficult to predict from the en¬ 
ergetics alone within Figure [3 ]d. One possibility is that 




Excited state order Excited state order Excited state order 


Figure 2: Wave function character density p (see text) for the excited states below the optical exciton in the eclipsed polyene 
cluster, for the (a) dimer, (b) trimer, and (c) tetramer, respectively. The different color schemes define the EXC (Scheme |^i), 
P + P _ (Scheme |]d), COV (Scheme I) and (TT) n (Scheme [Jl) contributions. Circled asterisks identify optical excitons with 
significant transition dipole couplings to the ground state (cf. Table [|. The lowest state in the tetramer is a CTX and not 
(TT) n . 


SF involves internal conversion from the 2 1 A ~-derived 
states to the CTX, followed by relaxation to the (TT) n . 
An alternate possibility involves branching of relaxation 
channels, whereby excitations decay to the CTX and the 
2 1 A~ in parallel, and only the former are involved in the 
SF. Lower yield of triplets is expected within the sec¬ 
ond scenario. Musser et al. have argued that in the 
polymer P3TV, in which the 2 X A ~ occurs below the op¬ 
tical 1 1 B+, SF occurs within the single-chain but directly 
without going through 2 X A~ state. EH The absence of 
involvement of the 2 1 A~ is in agreement with the pic¬ 
ture involving branching of the excitation decay process. 
The single chain mechanism however perhaps indicates 
bending and twisting of the long-chain polymer in the 
solution phase, such that different segments of the poly¬ 
mer are effectively decoupled. This can lead to CT be¬ 
tween decoupled segments. We postpone the discussion 
of the relatively slower SF in carotenoids compared to 
the acenes until later. 

Our calculations here place the (TT) n states below the 
COV 2 1 A~, which we believe is necessary for efficient bi- 
molecular SF. We have performed full Cl calculations for 
our Coulomb interaction parameters (U = 8.0 eV, n = 2) 
for octatetraene and decapentaene and have confirmed 
that E(2 1 A“)> 2xF?(Ti) in the single molecules. The 
relative locations of the 2 1 A~ and the (TT) n depend on 
both the correlation strength and the physical size of the 
molecule, as is easily understood from physical considera¬ 
tions. For strong Coulomb interactions, the single triplet 
exciton is strongly localized, and its physical “size” (cor¬ 
relation length, which is the expectation value of the sep¬ 
aration between the triplet-coupled spins) is small. In 
such a case, the 2 1 A“, which is a bound state of two 
triplet excitons, can be accommodated on a molecule that 
has length at least twice the size of the triplet exciton and 
E(2 X A“)< 2xi£(Ti). This is true for octatetraene and 
decapentaene for the standard Ohno parametrization of 


the PPP Hamiltonian (U = 11.26 eV, n = 1), as was 
found previously from exact single-molecule calculations. 
07108] We have confirmed these results also from full Cl 
calculations. For smaller and more realistic correlations, 
the size of the triplet exciton is larger, and in this case 
the state with two spin excitations has to be “squeezed” 
to fit into the short molecules; now confinement raises 
the energy, making £ 1 (2 1 A“) > 2 x F(Ti). A similar ar¬ 
gument was made previously for the bound state of two 
charge excitons (biexcitons) in linear chains, and numeri¬ 
cally demonstrated. m This effect can have an interest¬ 
ing consequence. On the one hand, with increasing chain 
length the energy of the lowest (TT) n state decreases rel¬ 
ative to the 1 X B+, making SF more exoergic and hence 
possibly more efficient. On the other hand, even with 
intramolecular Coulomb interactions, as the chain length 
is increased a length should be reached where E(2}A~) 
becomes less than 2 x E( Ti). Assuming that intermolec- 
ular Coulomb interactions have not raised its energy too 
much, from this length onward the lowest 2 X A“ and the 
lowest (TT) n state may be very close in energy, and this 
interference can be detrimental to SF since a smaller frac¬ 
tion of the excitations will end up as 1 (TT). In principle 
then, it is possible that SF efficiency reaches a peak at 
some intermediate molecular size. This would be inter¬ 
esting to examine experimentally. 

To summarize this subsection: (i) the optical state in 
polyene clusters is an optically allowed CTX within the 
eclipsed geometry but a Frenkel exciton in the mono¬ 
clinic crystal structure (several states are optically al¬ 
lowed with the slip-stacked geometry, but the state with 
the strongest oscillator strength is again a Frenkel ex¬ 
citon); (ii) the (TT) n states are the lowest excitations 
in the monoclinic structure and do not mix with the in¬ 
tramolecular 2e—2h 2 X A“, indicating that SF in polyenes 
is multichromophoric; (iii) neither do we find eigenstates 
that are superpositions of the EXC and (TT) n ; (iv) hence 
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4 chromophores (slip-stacked) 

■ COV 2e-2h sj 



8 10 12 14 16 18 20 22 24 26 

Excited state order 


Figure 3: Same as in Figure [2^i— c for clusters of four octate- 
traene molecules, with (a) slip-stacked and (b) monoclinic 
geometries. For the slip-stacked geometry, eigenstates do not 
have g and u symmetry anymore, and optically allowed states 
are strong superpositions of le—lh and 2e—2h excitations 
whose contributions have been shown separately. The excited- 
state spectrum for the monoclinic lattice predicts strong SF 
yield, in agreement with experiments. 


the SF in the monoclinic structure is mediated by the 
P + P - states, which are distinct eigenstates that occur 
between the EXC and the (TT) n ; and finally, (v) it is 
conceivable that SF yield in the polyenes peaks at some 
critical molecular size. We will see that the overall pic¬ 
ture changes somewhat in the acenes. 

Pentacene and Tetracene Crystals. Our calcula¬ 
tions here are for dimers (both tetracene and pentacene) 
and trimers (tetracene) of molecules occupying lattice 
sites of the characteristic herringbone crystal structure 
of acenes m shown in Figure |4^i. We chose intramolec¬ 
ular Coulomb and hopping parameters slightly different 
from the polyenes, in order to better reproduce known 
single-molecule results. Thus, while for the peripheral 
C—C bonds on the circumference of the molecules we 
chose the standard tC = —2.4 eV, for the internal bonds 
that are shared by the benzene nuclei we found better 
agreement with experiments by using = — 2.2 eV 
(the internal bonds are known to be weaker compared 
to the peripheral bonds EH). For the tetracene single 
molecule, QCI calculations with U = 5.0 eV and k = 2 


give the optical exciton and Ti at energies 2.23 and 1.21 
eV, compared to experimental values of 2.63 and 1.28 
eV in solution, |75] respectively (E( S*)< 2 x E( Ti) only 
in crystals). For the pentacene single molecule, the ex¬ 
perimental energies of the optical exciton and the triplet 
exciton are 2.12 and 0.78 eV, respectively. m The cor¬ 
responding quantities are 2.23 and 0.92 eV, respectively, 
from QCI calculations with U = 6.5 eV and k = 2 . In 
general, increasing U increases the energy of the singlet 
optical exciton but decreases the energy of the triplet 
exciton; the need to fit both the singlet and the triplet 
exciton energies requires compromising to intermediate 
Coulomb strength. These reduced Coulomb interaction 
parameters will therefore be used for our calculations re¬ 
ported below for the multiple-chromophore cluster. We 
find the 2 1 A~ at 2.57 (3.02) eV for the single molecule of 
pentacene (tetracene), i.e., higher in energy than the op¬ 
tical one-photon states by 0.34 (0.70) eV. Our results for 
the lowest two-photon state for pentacene contradicts the 
earlier work by Zimmerman et ah, who had found an in¬ 
tramolecular 2e—2h state below the optical state [33. and 
had suggested an intrachromophore SF. Recent experi¬ 
mental work on SF in amorphous films EH and solutions 
m of acenes both suggest the multichromophoric nature 
of SF. As before, our multiple-chromophore calculations 
are for k 1 - = k = 2 . We have performed these calcula¬ 
tions for a range of intermolecular tfj , with 0^|/3| < 0.2 
eV and appropriate relative angles between the p z or¬ 
bitals on neighboring molecules for the herringbone crys¬ 
tal structure in eq [7] (cf. Supporting Information for de¬ 
tails). 

Pentacene. We present our results first for pentacene 
since SF here is exoergic and the discussion is more 
straightforward. For pentacene, QCI can be implemented 
for a maximum of two chromophores only. In Figure^, 
c, we have given the same wave function analysis as for 
polyenes, for the dimers I—II and I—III of Figure [4^i, re¬ 
spectively, for (3 = — 0.2 eV. There is a significant differ¬ 
ence between the two dimers, although the lowest energy 
state in both cases is (TT)i. In the I—III dimer, be¬ 
cause of the very large distances between the molecules, 
CT plays no role whatsoever, and higher excitations are 
pure EXC and P + P - . This is different in the I—II dimer, 
where we see three different optically allowed excitations 
labeled S*, Sf and S 3 , respectively, in increasing order 
of energy. S* and S 3 are CTX in nature, with wavefunc- 
tion contributions from both EXC and P + P _ (S£ has 
minor P + P - contribution as well). Figure [4]T, e show 
the evolutions against (3 of all relevant wavefunctions of 
the I—II dimer as well as their energies and transition 
dipole couplings with the ground state. A complete and 
physical picture of the consequence of CT and its role in 
SF is obtained from the results presented in Figure [4]l , 
e, as we now discuss. At (3 = 0 , the dipole-allowed S 2 
and S 3 are degenerate, and S* is purely P + P - with no 
transition dipole coupling to the ground state. With in- 
























10 



Figure 4: (a) Herringbone lattice structure of acene crystals. Parameters R and 0 were taken from refs l32l and El respectively 
(the parameters are known to be very close for tetracene and pentacene, and hence the same set was chosen for both). QCI 
calculations are performed for the I—II and I—III dimers (both pentacene and tetracene) and for the I—II—III trimer for 
tetracene. The wave function analysis, as in Figure [2] and Figure 3] for the states proximate to the optical exciton in (b) for 
the I—II dimer and (c) for the I—III dimer of pentacene, respectively, (d) The relative contribution of P + P _ to the (TT)i, S^, 
S* and S* states of the I—II dimer in (b) as a function of /3 in eq [ t] (e) The magnitudes of the transition dipole couplings with 
the ground state of the allowed optical states for three different values of \/3\ = 0, 0.1, and 0.2 eV. Increasing intermolecular 
hopping brings the lowest optical exciton closer to the (TT)i and also increases its oscillator strength. 


creasing CT (|/3|), S* acquires EXC character and be¬ 
comes dipole-coupled to the ground state, while Sf loses 
oscillator strength as S 3 splits from it by acquiring P + P - 
contribution. Importantly, the (TT)i simultaneously ac¬ 
quires weak P + P - character, and S* also has moved 
closer to (TT)i. Moreover, the (TT)i state energy 1.79 
eV < 2 x E( Ti)= 1.84 eV, indicating the binding energy 
between two Ti triplets due to CT ((TT)i state energy 
is 1.83 eV at /3 = 0). Taken together, one therefore sees 
that in pentacene crystal there occurs an optical exciton 
that is very close in energy to (TT)i and that both the 
exciton and (TT)i have contributions from P + P - . We 
believe that the proximity in energies of S* and (TT)i, 
the shared partial P + P - character, and the exoergicity 
of the reaction lead to rapid generation of (TT)i and are 
behind the high SF yield in pentacene. An important 
difference from the polyenes therefore is that there are 
no intervening states here between the initial and final 
states. The experimentally observed Davydov splitting is 
0.14 eV. [76^ Theoretically, this is the splitting between 
S* and S 2 , and from Figure [4^ it seems (3 = —0.16 eV 
reproduces the experimental result closely. CTX charac¬ 
ter of the lowest le—lh excited states in pentacene has 
also been confirmed experimentally. EH [78] 

Tetracene. For tetracene, we were able to perforin QCI 
calculations for a cluster containing all three molecules I, 
II, and III in Figure |4^. However, as with pentacene, we 
first begin by looking at the I—II and I—III dimers whose 
energy-wave function analyses are presented in Figure [5^i, 
b, again for /3 = — 0.2 eV. Unlike pentacene, in tetracene 
the lowest X (TT) state is energetically higher than the 
optical exciton, and therefore it is taken to be the highest 
excited state in the energy window of our analyses. For 


the I—II dimer in Figure^, the (TT)i state has energy 
2.5 eV, very close to 2xF(Ti) = 2.42 eV. Figure^, b 
is qualitatively similar to Figure [4]d, c. Dimer I—III is 
again mostly irrelevant in the actual SF process, while 
as in pentacene in the I—II dimer the (TT)i state with 
partial P + P - character is proximate to an optical exciton 
(S 3 in this case) with significant P + P - character. Our 
calculated Davydov splittings (S*—S 2 and S*—S 3 gaps) 
this time are 0.15 and 0.22 eV, respectively, at /3 = — 0.2 
eV, compared to the experimental 0.08 eV, m indicating 
that the realistic \/3\ once again is slightly smaller. 

We show in Figure d the evolutions of the I—II 
tetracene dimer wave functions, energies, and transition 
dipole couplings to the ground states with /3. The results 
are similar to those for pentacene. For /3 = 0 optical exci¬ 
tation is to degenerate S 2 and S 3 excitons. With increas¬ 
ing /3, once again the S* loses P + P - character, acquires 
EXC character, and gains oscillator strength, while the 
dip ole-allowed S 3 blueshifts and acquires P + P _ charac¬ 
ter. Simultaneously the (TT)i state also acquires weak 
P + P - character. Thus, the mechanism of SF in tetracene 
is very similar to that in pentacene, in that direct tran¬ 
sition occurs from an energetically proximate optical ex¬ 
citon with significant P+P - contribution to (TT)\ with 
also partial P+P - character. One difference is that SF in 
tetracene crystals is probably still energetically uphill or 
at best nearly isoenergetic and, hence, less rapid than in 
pentacene. HJ m I n Figure^, we have given the wave 
function analysis for the I—II—III trimer. The new result 
is that now the lowest 1 (TT) state has a contribution 
from the (TT )2 component (triplets on I and III) that 
is more than three times the contribution from (TT)i 
(triplets on I and II), indicating a strong tendency to- 
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Excited state order |p| Energy (eV) Excited state order 


Figure 5: Wave function analysis for the states below the lowest (TT)i for (a) I—II and (b) I—III tetracene dimers, (c) Same 
as in Figure W- (d) Same as in Figure Et (e) Wave function analysis for the I—II—III tetracene trimer. (TT)i and (TT )2 
correspond to I—II and I—III dimers, respectively. 


wards triplet separation, in agreement with experimental 
observation. ng a similar tendency toward rapid triplet 
separation is also expected in the pentacene trimer. 

An interesting difference between pentacene and 
tetracene crystals is that in the former SF is from the low¬ 
est optically allowed singlet exciton Sf, while in the latter 
it is from the highest, S* in Figure pp. This has an in¬ 
teresting consequence for exciton dynamics in tetracene. 
Several research groups have recently pointed out that 
the “prompt” fluorescence in tetracene appears to orig¬ 
inate from two different emitters. p~6TFl9] Especially at 
lower temperatures a fluorescence channel opens up that 
appears to involve an emitter different from the one in¬ 
volved in SF. This lower energy emitter has variously 
been described as a “dark”, [16] “dull”, [17] or “trapped” 
m exciton. Our work suggests that this is simply the 
S* exciton and not a structural defect as suggested in 
ref fT6l Further experimental work is necessary to distin¬ 
guish between the two possibilities. 

Covalently Linked Acene Dimers. The high effi¬ 
ciency of SF in the acene crystals suggests similar high ef¬ 
ficiency also in covalently linked dimers of such molecules. 
The two triplets in this case are expected to occupy 
the acene fragments individually. Direct covalent link¬ 
ing leads to too strong CT, as in bianthryl, [80] which 
we have seen is detrimental to SF (see discussion of 
eclipsed clusters of polyenes above). Muller et al. there¬ 
fore investigated SF yield in three different bis (tetracene) 
molecules, linked through phenyl (tetracenes at para and 
meta positions) and biphenyl groups. [30, 31 In all cases 
the SF yield was insignificant, less than 3% for the best 
case. Muller et al. ascribed this to the redshift of the 
optical exciton in the bis (tetracene) molecule, relative to 
tetracene; since the X (TT) energy presumably stays the 
same in the dimers, the redshift of the optical exciton 
would imply larger activation energy to SF in this case. 
The very small shift in the optical exciton energy, as seen 
from the experimental solution absorption spectra (cf. 


Figure 1 of ref m suggests, however, that there is an 
additional and probably stronger reason for the tiny SF 
yield in the dimer (indeed, the authors have calculated 
the activation energies, and without assuming large dif¬ 
ferences in the prefactors in the Arrhenius rate equation, 
it is difficult to ascribe the gigantic difference in the SF 
yield to the difference in activation energy alone). 

We have theoretically investigated the bis (tetracene) 
molecule of Figure [ 6 ^ 1 , to understand the low SF yield in 
covalently linked dimers. Because of strong steric hin¬ 
drance, the phenyl group is expected to be rotated with 
respect to the tetracene segments. We assume that the 
twist angle 0 is the same in magnitude but in opposite di¬ 
rections for the two tetracene groups (cf. Figure^). The 
hopping integrals between the phenyl group and the pen- 
tacenes are taken to be tpp cos 6 , m where tpp = — 2.2 
eV is the standard hopping integral corresponding to a 
C—C single bond. We have performed QCI calculations 
for the bis (tetracene) molecule for varying 6 and show 
our results for the energies of the optical exciton and the 
X (TT) state of two triplets each located on a separate 
acene monomer in Figure^. Both energies rise rapidly 
with increasing 6. The very weak redshift in the experi¬ 
mental solution absorption energies of the bis (tetracene) 
molecules indicates that the realistic 0 in these molecules 
is large, perhaps even close to 90°, as is found also from 
ab initio calculations of the ground-state geometries. m 
For such large 0 , we have found that excitations are al¬ 
most entirely localized on the acene molecules without 
being delocalized also on the benzene linker. We have 
therefore labeled all intermolecular excitations with the 
suffix “AA”. The formal expressions for EXCaa, Pa^A’ 
and TaTa are the same as in eqs[5]and[6j with the only 
difference that the HOMO and LUMO now correspond to 
those of the acene units only. In Figure^, we have given 
the wave function analysis for the states below the low¬ 
est TaTa at 0 = 85°. The natures of the wave functions 
are now pure EXCaa, Pa^A’ an d TaTa as in Figure [ 5)3 
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for the I—III dimer of tetracene. We ascribe the low SF 
yield to this qualitative difference in the physical nature 
of the wave functions , rather than to the small difference 
in activation energy. Unlike in Figure [5^i, e, where an 
optically allowed CTX is proximate to the (TT)i with 
some P + P - component, states proximate to the TaTa 
are pure P + P - and are not optically accessible. Further, 
the P + P - states also separate the TaTa from the opti¬ 
cally allowed exciton. Figure^ also explains an appar¬ 
ently puzzling observation made by the authors of refl3T1 
The dimer according to the authors should have had two 
allowed exciton states, and the authors were perplexed 
that no such splitting of the absorption band was ob¬ 
served. As noted in Figure [6]:, the second exciton state 
is optically dark. 

It is interesting in this context to also extend our calcu¬ 
lations to the hypothetical bis (pentacene) molecule [36] 
since the issue of activation energy is not relevant in this 
case. We have performed the same calculations as for 
bis(tetracene) for bis (pent acene), assuming the same ge¬ 
ometry. In Figure [6jl, e, we have shown the ^-dependent 
energies and the wave function analysis for 0 = 85° for 
this case, respectively. Once again at this large 0 states 
are pure TaTa, EXCaa, and P^P^. Energetically, the 
two PtPX states now occur even above the two EXCaa 
states, only one of which is optically allowed (cf. Fig¬ 
ure^). The wave functions are similar to those for the 
I—III dimer of pent acene in Figure [4^, which we have 
noted contributes little to SF. We predict extremely low 
SF here too because of weak CT, even though the reac¬ 
tion would be exoergic. 

■ CONCLUSIONS 

In summary, we have performed QCI calculations of 
several different classes of molecules and their clusters, 
with the specific goal of arriving at a comprehensive phys¬ 
ical picture of SF in these systems. The attraction of the 
semiempirical approach is that it allows the inclusion of 
a very large number of basis functions, making it possible 
to treat all le—lh and 2e—2h excited states on equal foot¬ 
ing. This is necessary for arriving at the correct energy 
ordering of the excited states in the clusters, which we 
have argued is essential to understanding SF. The widely 
recognized inequality E( S*) > 2 x E( Ti) is a necessary 
but not sufficient condition for high SF efficiency. On the 
other hand, the PPP model QCI calculations allow us to 
explain the similarities as well as differences in SF be¬ 
tween three different classes of materials without making 
additional assumptions. 

For all the cases we have examined, we find that in- 
termolecular CT plays a strong role in SF. The precise 
role and extent of CT, are, however, strongly depen¬ 
dent on structure and morphology. In the acenes SF 
involves a one-step transition from a proximate singlet 


optical exciton to the 1 (TT) state, where both the initial 
and final states share P + P - character. The very high 
yields and rate of SF in pentacene are consequences of 
the unique electronic structures of both the optical ex¬ 
citon and the 1 (TT) states in this case. Our proposed 
mechanism here is very close to that given by Beljonne 
et ah, [38] for equilibrium separations between the pen¬ 
tacene molecules. We do not find a coherent superpo¬ 
sition of EXC, P + P - , and 1 (TT) with significant con¬ 
tributions from all three, and indeed such a superposi¬ 
tion is unexpected in the absence of direct coupling be¬ 
tween EXC and 1 (TT). Such a superposition has been 
proposed by Chan et al. [12] and has also been calculated 
by Beljonne et ah, for intermolecular separations smaller 
than equilibrium separations, when the large first order 
couplings between EXC-P+P - and 1 (TT)-P + P _ pairs 
can give large second-order coupling between EXC and 
1 (TT) if the P + P - state energy is not too large. Similar 
results were also found by Zeng et ah, [42] who had cal¬ 
culated the eigenstates of the pentacene I-II dimer in the 
reduced space of the six basis functions that are obtained 
when only the HOMO and LUMO of each molecule are 
retained in the calculation (the ground state with both 
LUMOs unoccupied, two EXC and P + P _ , and a single 
TT configuration as the basis components). Zeng et al.’s 
eigenstates are qualitatively similar to ours for equilib¬ 
rium geometry, except that they also find a state that is 
a strong superposition of EXC, P + P - , and 1 (TT) when 
the couplings between the molecules are made artificially 
large by either bringing them closer to one another or 
by rotating them. In neither of these cases the authors 
calculate the Davydov splittings for the increased inter¬ 
molecular couplings. It is our belief that the resultant 
Davydov splittings would be far larger than the exper¬ 
imentally observed quantities. This is substantiated by 
our determination that while with \f3\ ~ 0.16 eV in eq[7] 
we are able to reproduce the experimental Davydov split¬ 
ting of ~ 0.14 eV for pentacene even for a modest in¬ 
crease to |/31 = 0.2 eV in eq [7] calculated splittings exceed 
the observed values. We further believe that the mixing 
between EXC and X (TT) is reduced further for calcu¬ 
lations within the complete basis space (as opposed to 
the reduced space of six basis functions) since all three 
classes of eigenstates, EXC, P + P _ , and X (TT), now have 
contributions from higher energy excitations, CT involv¬ 
ing which is less efficient. Corresponding to each of our 
normalized many-body wave function calculated with a 
complete to nearly complete basis, it is possible to cal¬ 
culate the error in the computed wave function had only 
the six basis functions of Zeng et al. [42] been retained. 
This is simply the total sum of the contribution to the 
wave function by all other basis functions that are ig¬ 
nored within the reduced basis space. Proceeding in this 
manner we have determined that while the errors in the 
calculations of EXC and CTX in pentacene are 16% and 
12%, respectively, the corresponding error for 1 (TT) is 


13 



Figure 6: (a) Two tetracene molecules linked by a phenyl group. The acene molecules are rotated by angles ±0 with respect 
to the phenyl group, (b) The energies of the optical exciton and the lowest X (TT) (TaTa, see text) state as a function of 0 for 
tetracene. The horizontal dashed line is 2 x E(T±), where E( Ti) is the energy of the lowest triplet of isolated acene monomers, 
(c) Wave function analyses for the covalent dimer of (a), at 0=85° for tetracene. These wave functions are nearly the same for 
all 0 > 60°. (d) and (e) present equivalent analyses to (b) and (c), respectively, for pentacene. 


20 %. 

In contrast to the acenes, in polyenes with the mon- 
oclinic crystal structure, the transition is from a pure 
optical exciton to a pure 1 (TT) and is at least a two- 
step process, via P + P - states. Whether it is even a 
three-step process with the 2 X A“ as yet another inter¬ 
mediate or whether the relaxation of the singlet optical 
exciton occurs via two competing paths, one involving the 
2 X A“ and the other leading to SF, cannot be predicted 
from our calculations. In either case we expect slower 
generation rate of triplets through SF in the polyenes 
than in pentacene, in agreement with observations. El b 
|27|[52] Finally, errors incurred if the calculations of 1 (TT) 
states in polyenes are done in the reduced space of six 
basis functions are huge, close to 50%. This is because 
bonding (antibonding) MOs much lower (higher) than 
the HOMO (LUMO) make very significant contributions 
to the 1 (TT) here. 

The roles of CT are extreme in the polyenes with the 
hypothetical eclipsed geometry and the covalently linked 
acene dimers with realistic geometry. In the former, 
CT is too strong, and as a consequence the lowest ex¬ 
cited state is an excimer, rather than 1 (TT). Efficient 
SF is not expected here. There are many examples of 
real systems where excimer formation precludes SF. |80i 
The slip-stacked geometry is more favorable for SF than 
the eclipsed geometry because of weaker CT. In the co¬ 
valently linked acene dimers, in contrast, there is too 
little CT. The 1 (TT) states in bis (tetracene) are sepa¬ 
rated from the lower energy optical exciton by the opti¬ 
cally dark P + P - states at intermediate energy, preclud¬ 
ing transition from the optical exciton to the X (TT). This 
explains the very low SF yield here. We predict similar 
low SF yield also in the hypothetical bis (pentacene). 

The calculations we have reported are for a purely elec¬ 
tronic Hamiltonian, and vibronic couplings have not been 


included. The effects of including vibronic couplings can, 
however, be guessed based on the general observation 
that vibronic couplings lower the energies of molecular 
ionic states more than that for neutral molecules. The 
same principle applies to excited states, where vibronic 
couplings lower the energies of charge-separated ionic 
states more than charge-neutral covalent states. Further¬ 
more, the larger the electron—hole separation, the greater 
the energy lowering. Thus, energies of polaron pairs and 
charge-transfer exciton states, with the charges on sepa¬ 
rate molecules, will be lowered more than those of the in¬ 
tramolecular exciton, and the latter’s energy will be low¬ 
ered more than triplet—triplet states. From Figure 
this would imply even larger SF rate for pentacene where 
S* becomes more proximate to X (TT). Conversely, S* in 
Figure^ becomes less proximate to the X (TT) above it 
in tetracene, and the SF rate is reduced slightly. In the 
monoclinic carotenoids, given the very large numbers of 
states involved not much change is expected. In the case 
of the bis (tetracene), vibronic couplings will not make 
any difference, as the states closest to the triplet—triplet 
are pure polaron pairs that are not optically accessible. 

Finally, our work here has focused on the mecha¬ 
nism of SF alone and not whether the molecules we 
have investigated will be particularly useful in enhanc¬ 
ing photovoltaic performance. A necessary condition 
for the triplets undergoing charge dissociation at a 
donor—acceptor interface is that the donor—acceptor ex- 
ciplex has even lower energy. [3] This condition is cur¬ 
rently satisfied by relatively few heterostructures, the 
most well-known example of which is the pentacene—Ceo 
solar cell. The low triplet energy in this case is the 
reason behind the very small open-circuit voltage, 0 - 
E] which reduces the performance of the solar cell even 
though the number of charge carriers is large. Organic 
donor molecules with significantly larger triplet energy 
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than acenes are therefore of interest, and we are aware 
of experimental investigations of polycyclic aromatic hy¬ 
drocarbons (PAHs) with large optical and triplet gaps as 
candidates for efficient SF. [82 We have recently per¬ 
formed theoretical work on excited state orderings in 
large PAH molecules. [83] The approach we have de¬ 
veloped here can in principle be extended to clusters of 
such PAH molecules. This is the goal of future work. 
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■ NOTE ADDED IN PROOF 

After this manuscript was accepted for publication, we 
were notified about the new publication by T. C. Berkel¬ 
bach et ah, [84] in which the authors have also discussed 
CT character of the low-lying excited states in crystalline 
pentacene. 
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1. Details of QCI calculations 

Number of C atoms, active MOs and dimensions of the final QCI matrices for octatetraene (cf. Figure la), 
decapentaene, pentacene, tetracene (cf. Figure 4a) and covalently linked dimers (cov-dim) of tetracene (cf. Figure 6a) 
and pentacene multiple chromophores are all given in the following table: 




octatetraene: 


decapentaene: 

pentacene: 

tetracene: 

tetracene (pentacene) cov-dim: 

Chromophore number: 

2 

3 

4 

2 

3 

2 

2 

3 

3 

C atom number: 

16 

24 

32 

20 

30 

44 

36 

54 

42 (50) 

Active MO number: 

16 

18 

16 

20 

18 

20 

20 

18 

20 

Matrix size: 

294185 

772435 

294185 

1833276 

772435 

1833276 

1833276 

772435 

1833276 


Table SI: Details of QCI calculations for multiple chromophores of polyenes and acenes studied in the main article. 


2. Monoclinic lattice structure and intermolecular hopping integrals for four oc¬ 
tatetraene chromophores 

According to refs [SI, S 2 ], diphenylpolyenes usually crystallize in monoclinic and orthorhombic lattice structures. 
The monoclinic structure was adopted for this work in which more efficient SF compared to the orthorhombic 
structure was reported. [S3] We take an octatetraene chromophore originally lying in the xy plane along the x axis 
with its center on the origin. We rotate this chromophore around the x axis by 7 = ±30° followed by another rotation 
around the 2 axis by 0 = 42.30°, resulting in two new chromophores. [SI] A copy of the one rotated by 7 = +30° is 
translated by distance a = 6.23A along the x axis while the other rotated by 7 = —30° and a copy of it by distances 
\a along the x and along the z axes with b = 7.45A (cf. Figure lc). [S2] Using eq 6 in which /3 = — 0.2 eV 
and h m i n = 3.5A as chosen for the eclipsed geometry (cf. Figure la) and also considering the angles between the 
p z orbitals and the lines connecting them on two separate chromophores, the largest calculated t 1 - value was —0.04 eV. 

3. Intermolecular hopping integrals for acene chromophores in the herringbone lattice 
structure 

For acene chromophores possessing herringbone lattice structure with parameters given in Figure 4a, we employed 
eq 6 with /3 = — 0.2 to —0.08 eV in steps of 0.02 eV (cf. Figure 4d, 5c) and h m i n = 3.50A to compute t 1 - values. 
Taking account of the relative angles between the p z orbitals, only two intermolecular hopping integrals, i.e., t± and 
£ 2 “, become greater than 0.07 eV, our intermolecular hopping integral threshold, for the largest /? = — 0.2 eV as shown 
in Figure S 2 for two pentacene chromophores. Our chosen range of (3 is consistent with the weak to intermediate 
CT effects between I—II and II—III chromophores in Figure 4a in comparison with experiments. [S4, S5] Within 


SI 




Figure SI: (a) Side and (b) top views of intermolecular hopping integrals among four octatetraene chromophores in the 
monoclinic lattice structure with lattice parameters given in Figure lc. Angles 61-4 represent the relative orientations of the 
p z orbitals on two separate chromophores. 


this range, our calculated t 1 - values as presented below between the sites numbered for two pentacene chromophores 
depicted in Figure S2: 

P = -0.2 eV: 

ti = -0.13 eV: (24,5), (30,9), (34,13), (38,17), (42,21). 

= -0.07 eV: (24,4), (24,6) (30,6), (30,10) (34,10), (34,14) (38,14), (38,18) (42,18), (42,22). 


and 

P = -0.18 eV: 

ti = -0.12 eV, t£ = -0.07 eV 

P = -0.16 eV: 

tj- = -0.11 eV, w 0 

P = -0.14 eV: 

tf = -0.09 eV, « 0 

P = -0.12 eV: 

ti = -0.08 eV, w 0 


S2 






















Figure S2: View of intermolecular hopping integrals between two pentacene chromophores in the herringbone lattice structure 
with lattice parameters given in Figure 4a. Angles 61-4 represent the relative orientations of the p z orbitals on two separate 
chromophores. 


P = -0.1 eV: 

ti = -0.07 eV, « 0 

P = -0.08 eV: 

ti = -0.06 eV, f 2 x «0 

where |t^| < 0.07 eV ~ 0. The same t 1 - applies to the equivalent sites of two tetracene chromophores. 

4. 2 X A“ state energy in the eclipsed, slip-stacked and monoclinic arrangements 

of hexatriene chromophores 

QCI results for the 2 X A“ state energy of hexatriene chromophores are presented below as they are added one 
by one in the eclipsed geometry (cf. Figure la) and also when four chromophores are arranged together in the 
slip-stacked and monoclinic geometries (cf. Figure lb, c). The number of MOs included in QCI calculations is also 
specified (energies are in eV): 


Number of chromophores: 1 

2 

3 

4 

4 

4 

Number of MOs included in QCI: 6 

(eclipsed) 12 

(eclipsed) 18 

(eclipsed) 16 

(slip-stacked) 16 

(monoclinic) 16 

Ei = 3.98 

E± = 4.13 

E\ = 4.45 

Ei = 4.76 

Ei = 4.67 

Ei = 4.93 


E 2 = 4.16 

E 2 = 4.51 

E 2 = 4.81 

E 2 = 4.70 

E 2 = 4.93 



E 3 = 4.52 

E 3 = 4.87 

E 3 = 5.12 

E 3 = 4.93 




E 4 = 4.98 

E 4 = 5.16 

E 4 = 5.04 


Table S2: QCI results for the energy of 2 X A g states in the eclipsed, slip-stacked and monoclinic arrangements of hexatriene 
chromophores. 


As it is readily seen, 2 1 A g state energy increases with the number of chromophores. 
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